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Abstract. We consider principal moments of inertia of 
axisymmetric, magnetically deformed stars in the context 
of general relativity. The general expression for the mo- 
ment of inertia with respect to the symmetric axis is ob- 
tained. The numerical estimates are derived for several 
polytropic stellar models. We find that the values of the 
principal moments of inertia are modified by a factor of 2 
at most from Newtonian estimates. 
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1. Introduction 

Various features of pulsars have been disclosed with grow- 
ing observational data. Most pulsars have stable pulse 
shapes, and spin down steadily with typical time-scales of 
order several Myr (see e.g. Meszaros 1992j ). However, devi- 
ations from the linear spin-down trends have also been ob- 
served on shorter time-scales in some pulsars. These devia- 
tions convey information about the surroundings, internal 
structure and dynamics of neutron stars. Fo r exa mple, the 
pr esence of glitches (see e.g. Shemar & Lyne 1996 ; Wang et 
al. 2000) indicates the internal structure that two or more 
poorly coupled components coexist in the stars, while tim- 
ing noise (see e.g. D'Alessandro 1995) is believed to be 
due to random movements of the fluid in pulsars. Besides 
these short time-scale instabilities, the spin-down fea- 
tures indicating the precession due to stellar deformation 
have also been reported. Two anomalous X-ray pulsars 
(AXPs), 1E1048. 1-5937 and 1E2259+586, were obs erved 
to indicate irregular spin-down (see e.g. Mereghetti 1995; 
Bayk al fc S wank |1996[ Baykal et al. |1998| ; Oosterbroek 
et al 



1998). AXPs are a subclass of X-ray pulsars with 



period in a very narrow range (6-12s) and period deriva- 
tives in a range (10~ 12 — 10 _11 ss _1 ) (see e.g. Mereghetti & 
Stella |l995 ), and those are most likely candidates of mag- 
netars (Duncan & Thompson 1992| ; Thompson & Dun- 



1993; 1995 



1996) along with s oft ga mma repeaters 
(S GRs) (see e.g. Kouveliotou et al. 1998 Kouveliotou et 
al. |1999| ). As discussed by Melatos p999| ) , the wobbles 



in their spin-down rates may be interp reted by an ef- 
fect called radiative precession (Melatos 200C), which is 
related with an oscillating component of electromagnetic 
torque . Ano ther precessing pulsar PSR B1828-11 (Stairs 
et al. 2000 ) was reported very recently. This gives the 
very clear observation of free precession of an isolated 
neutron star. These deformed, precessing objects are ana- 
lyzed by solving the Euler equation of motion in the form 



■ SijkP fi ' fij = Ni, where is the inertia ten- 



sor, fl l is the angular velocity and Ni is the torque acting 
on the object. When we take into account the electro- 
magnetic torque by a rotating magnetic dipole (Davis & 
Goldstein 1970 ; Goldreich 1970 ), the above-mentioned ra- 
diative precession can be found for magnetically deformed 
stars (Melatos [1999; 



200C). Furthermore, since magnet- 



ically deformed, rotating stars emit gravitational waves, 
we can consider the gravitational radiation reaction torque 
(Bertotti & Anile |1973| ; Cutler & Jones |200C| ). The grav- 
itational backreaction damps the wobbles on a time-scale 
proportional to [(I z — I x ) /I x ]~ 1 f° r wobbling, axisym- 
metric rigid bodies, where Ii is the principal moments of 
inertia. Thus, the moments of inertia play a significant 
role in the analyses of pulsar precession. The estimates of 
the moments of inertia have been done in the context of 
Newtonian gravity so far. However, neutron stars are fully 
general relativistic objects, and it is important to take into 
account general relativistic effects. Therefore, we now dis- 
cuss the principal moments of inertia of magnetically de- 
formed stars in the context of general relativity. 

In Newtonian gravity, once we have the the mass dis- 
tribution of an object p(x), the moment of inertia with re- 
spect to any axis can be calculated from I = J p(x)x 2 d 3 x, 
where \ denotes the length from the axis. In general rel- 
ativity, only for axisymmetric objects, the moment of in- 
ertia with respect to the symmetric axis can be well de- 
fined. For this definition, we need slow rotation of the 
objects. The slow rotation ensures that the angular mo- 
mentum J is linearly related to the angular velocity f2, 
i.e. J = Ifl. Here, I defines the general relativistic ver- 
sion of the moment of inertia. The moments of inertia of 
relativis tic, sp herically symmetric stars were discus sed by 
Hartle (|l967| ) and Chandrasekhar & Millar (|l974| ). The 
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approximate expression valid for various realistic equa- 
tions of state was also derived by Ravenhall & Pethick 

P9D - 

In this paper, we discuss the principal moments of iner- 
tia of magnetically deformed stars which are axisymmetric 
with respect to the magnetic axis. Deformation of rela- 
tivistic magnetized stars was studied both in the numeri- 
cal approach (Bonazzola et al. |1993| ; Bocquet et al. 1995) 
and in the analytic approach (Konno et al. 1999). Based 



on the la tter, analytic approach developed by Konno et 
al. ( 1999| ), we proceed with t he discussions. In the previ- 
ous work (Konno et al. 1999| ) , stellar magnetic fields were 
regarded as corrections to non-rotating, spherical stars, 
and the magnetic deformation of stars was formulated us- 
ing a perturbative approach. As mentioned above, we now 
take into account the slow rotation of the deformed stars 
on the symmetric axis to define the moment of inertia. 
The formulation of this configuration is given in §||. The 
general expression for the principal moment of inertia with 
respect to the symmetric axis is also derived. The numeri- 
cal estimates are obtained for several stellar models in §0. 
In §^|, we discuss the other principal moments of inertia. 
Finally, we give conclusion in §||. Throughout this paper, 
we use units in which c = G = 1 . 

2. Magnetically deformed rotating stars 

We consider a slowly rotating star which is subject to 
quadrupolc deformation du e to a dipolc magnetic field. As 
described by Konno et al. ( 1999 ), magnetic fields can be 
treated as corrections to a background, spherically sym- 
metric star, i.e. B ~ O(eb)- We now assume that the star 
slowly rotates on the magnetic axis with a uniform angular 
velocity ~ 0(en), which is also regarded as a pertur- 
bation. In this paper, we take into account the rotational 
corrections up to first order in £q. The metric describing 
such a star can be given by 

s"W [i + 2 (fco(r) + h 2 (r)P 2 (cos6))] dt 2 

2e A M " - 
(mo(r) + m2(r)P 2 (cos0)) 



ds 2 = 



1 



dr z 



r 2 {1 + 2k 2 (r)P 2 (cos9)) d9 2 
(l + 2fc 2 (r)P 2 (cos60) 

W 3 (r) dP 3 (cos9) 



? • 2 

-r sin 



d(j>- VFi(r) 



sin 9 



d0 



n 2 



dt 



(1) 



where Pi is the Legendre polynomial of order I, and v 
and A describe a background star. The corrections ho, h 2 , 
mo, m 2 and k 2 of second order in eb correspond to de- 
viation from the spherical shape. If we replace Wx and 
W3 with zero in Eq. ([!]), this expression reduces to the 
metric already investigated by Konno et al. (1999). The 



newly appeared functions Wx and W3 lead to frame drag- 
ging due to the magnetically deformed, rotating star. In 
order to involve the effect of deformation, let these func- 
tions include corrections up to order £qE b . Showing order 



explicitly, we can write down W± and W 3 in the forms 



Wx = w + W, 



(2) 



W 3 = W: 



(2) 



(2) 
(3) 



where u ~ O(en) and (W x (2) , W^) ~ 0(e n e B ). The form 
of the terms including W± and W3 in Eq. (|]) corresponds 
to quadrupole deformation of the star. These functions are 
very analogous to the rotational corrections up to third 



order in eq (Quintana 1976; Kojima & Hosonuma 1999). 

The stress-energy tensor of the star which is composed 
of the perfect fluid endowed with a dipole magnetic field 
has the form 



1 

4tt V" 4 



• — ( /- /,X /va - -F^F^^i 



(4) 



where the density p and the pressure p are expanded as 



p = p {r) + (p 2 a{r) + p 22 {r)P 2 ) , 
P = Po(r) + (p 2 o(r) +p 22 {r)P 2 ) . 



(5) 
(6) 



Here, p 2 o, p 22 , p 2 Q and p 22 are second-order quantities in 
Eb- The non- vanishing components of the 4- velocity w M 
are given by 



-* [1 - {ho + h 2 P 2 )} 



= VLuK 



(7) 



Furthermore, the Faraday tensor F^ v can be giv en by 
the 4-potential A u writ ten as (Muslimov & Tsygan 1986 ; 
Konno & Kojima |2000|) 



A-n = a t o(r) + a t2 (r)P 2 ,Q, 0, a^r) sin 



de 



(8) 



where the 0-component corresponds to a dipole magnetic 
field, and the i-component denotes the electric field in- 
duced by stellar rotation. Therefore, we have ~ 0(eb) 
and (a t o,a t2 ) ~ 0(£n£s). 

The above-mentioned functions except W\ and W3 
were investigated i n det ail by Konno et al. ( 1999 ) and 
Konno & Kojima (200C) (see Appendix [a| for the brief 
summary). The differential equations which W\ and W3 
obey can be obtained from the (£0)-component of the Ein- 
stein equation, 



LA. 

r 4 dr 

l_d_ 

r 4 dr 



4 .dWi 
r J—\ — 
dr 



j_do_ 

r dr 



-j s 



S2 
5 



(9) 



.dW 3 



r J- 



dr 



4dj 10 .S 2 

je * W 3 =J — , (10) 

r dr r z I 5 



where j and W \ arc defined, respectively, as 



J 



Wx = n~Wx = w- w[ 2) {w = n - u) 



(ii) 



(12) 
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These differential equations are solved order by order. The 
differential equation for zu of order eq, in which both 



Si and 5*2 vanish, was solved by Hartle (1967). We are 
now interested in the differential equations of order EnS B . 
The functions Si and S% appearing in the source terms in 
Eqs. (|]) and JlC|), in general, include W\ and W3. How- 
ever, when we consider S\ and S2 up to the order of our 
interest, these functions include ft and 10 only as shown 
below. The source terms of order ens B are given by 



Si 



-h Q m 

r 



A') wmo 



■16-Ke x w (p 2 o +P20) + 



16e } 
Ik 4 



3r 2 



co - 



3r 2 



8e^ 
r 



-a<pa t 2 - -^a^a t0 , 



(13) 



-m 2 



32e 



r e A (i/' + A') mrri2 



■P22) ■ 
16e A 



3r 4 



-a0a t2 2 a <P a t2' 



(14) 



where the prime here denotes differentiation with respect 
to 7'. These functions can be calculated by using the results 
for the density and pressure (^20, ^22,^20,^22), the metric 
functions (ho, h%, mo, mi, te) and the potential functions 
(q<j), Qto, a ti), which were already calcu lated by Konno et 
al. ( |1999|) and Konno & Kojima (|2000| ). 

The differential equations ([))) and (|l(]) can be solved 
numerically by imposing boundary and junction condi- 
tions. The boundary conditions are summarized as 



W\ — ► const, W 3 — > 



as 



Furthermore, we i mpos e the junction condi tion g iven b y 
(O'Brien & Synge |1952| see also Eqs. ( |A.10| ) and flA.llD ) 



W'l^fm = W'\ 



(17) 
(18) 



-g(0) ) 

1+^(0) = vy a \-£(o) , 

where a takes 1 or 3, and denotes the surface of the 
background star. 

Before deriving numerical solutions, it is worthwhile to 
investigate the behavior of W\ and W3 at large r in detail 
using Eqs. (g), @, (|l|) and (Q), because this inspection 
gives the expressions for the angular momentum and the 
moment of inertia of the star. 

At large r, we have 



A. 



0. 



i.e. 



From Eqs. 
iors 

Si « ^ 



inl- 
and (0) 

1 



S 2 oc 



(19) 

we derive the asymptotic behav- 



(20) 



Using these expressions, from Eqs. (||), ( |l0| ) and (|l6|), we 
obtain 



wf oc 



TF- 



(2) 



Hence, we can put 

Wi^-^-<\ 



(21) 



(22) 



where J is the angular momentum of the deformed star, 
and w/p ~ 0(l/r 6 ) is a function of order eoe^. We now 
integrate Eq. (g) after multiplying both sides by r 4 . Us- 
ing the above expression (|2^), we can obtain the general 
expression for the principal moment of inertia 



r A — —dr 
dr ft 



+ - 



ft 



1 

6ft 



dr 

R 

4 • 

r J 





dr 



Si- 



dr, 



(23) 



where i? is the radius of the background star. The moment 
of inertia of the background, spherically symmetric star 



j(o) i s gi ven by (Hartle |g67| Ravenhall & Pethick |1994[ ) 
2 f R ^ cfl? 



j(o) 



r "d?n dr - 



3 jo 

Therefore, we obtain 



(24) 




r(2) 



3*^ , 

r — -= — dr 



dr ft 



(25) 



where we have used the decomposition I z = +1^ . The 

magnetic modification of the principal moment of inertia 

(2) 

Iz is the second order quantity in e B - The numerical 



estimates of I 
next section 



(2) 



for several stellar models are given in the 



3. Numerical estimates of the moment of inertia 

Using the boundary and junction conditions mentioned in 
the last section, we can derive the numerical solutions for 
W\ and W3 for any stellar models. The numerical results 
for the magnetic correction of the principal moment of in- 
ertia /i 2 ' can also be obtained by using Eq. ( |25| ) . We now 

(2) — 
show the results of I z obtained for several polytropic 

stellar models. In these calculations, we adopted a differ- 
ent condition from that of the previous study (Konno et 
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Fig. 1. The magnetic correction of the principal moment Fig. 2. Comparison between A/i and A/2 for n = 1. A/i 

of inertia if 1 plotted as a function of M /R. The values are and AI 2 are normalized by / M and plotted as a function 

normalized by ^t 2 /M, and the polytropic index is denoted of M/R. 
by n. 



al. 1999), in which sequences with constant central den- 
sity were investigated. In the current case, in order to clar- 
ify the correspondence between non-magnetized and mag- 
netized stars having same mass, we use the condition in 
which the total mass of the star does not change through 
the perturbative approach. This is also because the mo- 
ment of inertia can be modified significantly by the mass 
shift rather than the magnetic deformation. In our for- 
mulation, this condition is accomplished by imposing the 
boundary condition that mo vanishes at infinity. 

Figure [l] displays the magnetic correction 1^ as a 
function of the general-relativistic factor M/R. The values 
are normalized by the typical value fi 2 /M . As a simple ex- 
ample, we now discuss quadrupole deformation of a fluid 
body in the Newtonian limit. If the deformed body has 
constant density, which corresponds to n = 0, then we can 
derive the result 1^ = 10/i 2 /(3A /) fro m the estimat e of e l- 
lipticity 25n 2 /(2M 2 R 2 ) (Ferraro |l954i Konno et al. |2000D . 
However, as seen in Fig. our numerical result for n = 
is different from this simple estimate. This is because the 
perturbed star does not have constant effective density 
by the added perturbation, even though we assume the 
background star with constant density. This reason can 
be understood by seeing the differential equation for mo, 
i.e. Eq. (A3). The derivative of mo is related to the effec- 
tive density including electromagnetic energy. In the case 
of n = 0, although the first term on the right-hand side 



in Eq. (A. 3) vanishes, the remaining terms do not vanish 
and are non-trivial functions. It follows that the effective 
density is not a constant. Thus, our results include the in- 
ertia of electromagnetic fields as well as mass. Therefore, 
our result shown in Fig. fj cannot be compared with the 
above simple estimate. 

Concentrating on the general relativistic effects on the 
magnetic correction 1^ , we can find from Fig. [l] that the 



(2) 

values of I z for each stellar model become large with the 
general relativistic factor M/R. The increments are 50% 
at most. However, the rates of increase may be neglected 
except for the case of n = 1.5. 



4. The other components of moments of inertia 

Next, we discuss the other principal moments of inertia, 



i.e. I x and ly' . In the context of general relativity, the 
definition of these quantities is associated with a difficulty 
in the concept. As seen from Eqs. ( ^2|) and the notion 
of moments of in ertia i s related with quantities at infinity 
(see also Geroch 1970 Hansen 1974 ). If we consider the 
rotation of the star on the x-axis or the y-axis, it pro- 
duces the radiation of electromagnetic and gravitational 
fields. Thus, the exterior space-time is not stationary, but 
radiative. There is no rigorous way to define the moments 
of inertia in such a radiative space-time. The concept of 



r(2) 



li 2 ' and Iif' itself may be meaningless to a considerable 
extent. 

Therefore, only approximate expressions are available. 
It would be useful to extend Newtonian expressions, since 
it seems that most of Newtonian features still survive for 
neutron stars specified by the general relativistic factor 
M/R ~ 0.2. We cannot prove the validity mathematically, 
but will test some empirical relations. For this purpose, let 
us recall relations using the principal moments of inertia 
and the ellipticity which hold for an incompressible fluid 
in the context of Newtonian gravity. These are helpful in 
seeking the other principal moments of inertia. First, if the 
fluid body is subject to quadrupole deformation, then the 
relation between the principal moments of inertia exists 



r(2 



7(2) 



r(2) 
v 



.1/(2). 



(26) 
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Using this relation, the difference between the principal 

(2) (2) 

moments of inertia I z and I x is given by 



A/ 1 = /i 2 )-4 2 ) = |/i 2 ). 



(27) 



In the different way, this difference can be expressed also 
by using the ellipticity, which is defined by the difference 
between the equatorial radius and the polar radiu s (see 
e.g. Chandrasckhar & Miller |l974j Konno et al. |f999 ; 
Konno et al j2000| ), 

x (ellipticity). (28) 



Ah 



- r(o) 



Thus, the difference is derived by the two methods. There 
is no reason that the two expressions coincide, but empir- 
ically the agreement is good. Note that these expressions 
( p7j ) and (|||) coincide exactly in this case. 

If we assume that the relation (^) is applicable to gen- 
eral relativistic cases, it seems that we should simply apply 
the relation (pq) to the results shown in Fig. |] and derive 

(2) (2) 

I x and Iy . However, in order to utilize this relation, the 
magnetic corrections must be purely quadrupole contribu- 
tions. As seen from the metric (Q), the magnetic deforma- 
tion includes monopole parts as well as quadrupole parts. 
Thus, our results in Fig. ^ also include monopole parts. 
Hence, we have to subtract monopole contributions from 
Iz 2 ^ to derive I x 2 ^ and I y 2 \ It is, nevertheless, not so easy 
to extract monopole contributions accurately from these 
magnetically deformed stars. 

However, there is no guarantee that all the results in 
Fig. ^] include monopole contributions dominantly. If there 
is the case in which the monopole part can be neglected, 



then we can estimate I x 2 ^ and Iif' in the case, by assum- 
ing that the relation ( ptf ) holds approximately. That case 
would also provide some extrapolation for general rela- 
tivistic effects on I x 2 ^ and Iy in the other cases. In order 
to seek such a case, we now compare A/i derived by sim- 
ply using the results in Fig. [l] with A/ 2 calculated from the 
ellipticity, which was already estimated for several stellar 
models (Konno et al. 2000| ) . In the case that Aii is al- 
most consistent with AI2 about the values and tendency, 
it seems that the monopole part can be neglected, and 
we may use the relation (p6|). From the comparison, we 
find that A/i is almost consistent with A/2 in the case of 
n = 1. Figure |2| displays the comparison between A/i and 
A/2 in this stellar model. The two curves coincide within 



r(2) 



10%. Hence, we can estimate I x 2 ^ and Iy*' using Eq. (|26|) 

(2) (2) — 

in this stellar model. Since the values of I x and I y are 

(2) 

simply given by multiplying I z by a factor of —1/2, the 



r(2) 



(2) 1 

changes of the absolute values of I x and Iy* 1 due to gen- 
eral relativistic effects are specified by the same factor as 
in the case of 1^ (see Fig. ||). Consequently, we derive 
very similar result to that of I z 2 \ However, I x 2 ^ and I y 2 ^ 

decrease with the general relativistic factor M/R, while 
f 2) 

Iz increases. We expect that similar features exist also 
in the other stellar models. 



r(2) 
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Fig. 3. The other components of the principal moments 
of inertia I x 2 ^ and I y 2 \ which are derived for n = 1. 



5. Conclusions 

Recent observations indicate precession of some pulsars. 
In order to understand the dynamics of these pulsars in 
detail, we have to know the moments of inertia, which play 
a crucial role in the Euler equation of motion. Motivated 
by this consideration, we have discussed the general rel- 
ativistic effect on the moments of inertia of magnetically 
deformed stars. By developing the formulation by Konno 
et al. ( |1999D , we have considered the slow rotation of mag- 
netically deformed stars on the symmetric axis in order to 

define the moment of inertia. The general expression for 

(2) 

the magnetic correction I z to the moments of inertia of 
spherically symmetric stars was obtained, and the numer- 
ical estimates for various polytropic stellar models were 
also obtained. The other components of the moments of 

(21 (2) 

inertia I x and I y were discussed from the extension 



of relations in Newtonian gravity. From those results, we 
found that each principal moment of inertia is modified by 
a factor of 2 at most due to the general relativistic effect. 
Nevertheless, it seems that the general relativistic effect 
does not affect precessing motion noticeably in the case of 
free precession of pulsars. However, it is not clear whether 
or not the electromagnetic radiation reaction torque act- 
ing on the pulsars change due to the general relativistic 
effect. If this torque is modified by the general relativistic 
effect, then the period of the wobbles would also be mod- 
ified. Therefore, this point should be made clear in future 
investigation. 
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Appendix A: Summary of the perturbative 
approach 

we now review the perturbative approach developed by 
Konno et al. ( |l999| ) and Konno & Kojima ( [2000D 



A.l. Equations of order Eb 

First, we consider the magnetic field described by a^, 
which has the lowest order. This function obeys the fol- 
lowing equation derived from the Maxwell equation, 

+ i (✓ - V) e" ^ - = -A, U . (A.1) 

The electric cu rrent must satisfy t he intc grability con- 
dition (Ferraro 1954 ; Bonazzola et al. 1993 ) 



U = c r 2 (po+Po) 



(A.2) 



where Co is a constant of order eb- This corresponds to 
the c ondi tion in which the star is in equilibrium. Equa- 
tion ( A.l ) should be solved by imposing the condition in 
which a0 vanishes at the stellar center and infinity and is 
connected smoothly at the stellar surface. The exterior so- 
lution can be obtained analytically (Ginzburg & Ozernoi 



1965; Petterson 1974), while the interior solution was ob- 



tained numerically (Konno et al. 1999 ). 
A.2. Equations of order e 2 b 

Next, we consider the effect of magnetic stress, which is of 
second order in eb ■ This effect arises from the electromag- 
netic part of the stress-energy tensor (^) . The second-order 
metric functions ho, hi, mo, m 2 and k 2 can be obtained by 
solving the following two sets of differential equations and 
one algebraic equation derived from the Einstein equation 
G^ v = 87rT^, and the equation of motion Tf£,.„ = 0, 



drrio 
dr 



dho 
dr 



dv 2 
dr 

dh 2 
dr 



-47rr 2 ^(p +Po) (ho -a) 
Po 



■\ e ~ x K>) 2 



2 2 



8tt P 'o . 

Y7o a ^' 



(A.3) 



1 

i 

3~7 



H e rn - Airre (p + Po) (ho - ci) 



,A„2 



Sn 



3r3 C a ^ + ^ e " c W' 



(A.4) 



-v'h 2 + \e x (J- + V -\ (a;) 2 + -^a^, (A.5) 



;V 2 



e\ . 2 1 - e> 

8tT — ( P o+Po) + 2 



3rV 



. rv . , 

w v 1 + — 



+ Q-- A + ^) K) 2 + ^a^, (A.6) 



TO 2 



(A.7) 



where C\ is a constant fixed by the junction condition at 
the stellar surface, and v 2 is defined as v 2 = h 2 + k 2 . 
The differential equations ( |A.3| )-( [A~6| ) should be solved 
by imposing the boundary conditions 

mo, h 2 , k 2 — + 0, ho — > const as r — ► 0, (A. 8) 



^Oj ^2 - * 0, toq — ► const, 



1 



— (a > 3) 



as r — > oo. (A. 9) 

Here, note that the value of too at infinity corresponds to 
mass shift associated with the perturbation. Furthermore, 
the junction condi tion a t the stellar surface is written as 
(O'Brien & Synge |l952| ) 

ff/xwl+j = 9iiv\-£ (P>>" = t>r,0,<f>) , (A.10) 
9\o,t\+£ = 9\oA-£ (\v = t,0,<p), (A.ll) 
where £ denotes the boundary surface of the star. 

Once we have the above metric functions, we can derive 
the second-order corrections of pressure and density from 
the relations 

2 

Po)h + —HH + c i (Po 



P20 



P22 



P20 



P22 



- (Po 

- (Po 
Po 

—P20, 

Po 
Po 

—P22- 

Po 



Po) h 2 



3r 2 
2 

3^2 HH* 



Po) , 



(A.12) 
(A.13) 
(A.14) 

(A.15) 



Therefore, the stellar 
magnetic deformation 
Eqs. 



structure which is subject to 
can be determined by solving 
S)— (A.6). The exterior solution for the metric 
functions was obtained analytically, while the interior so- 



lution was calculated numerically by Konno et al. (1999). 
The fact that the evaluation of ellipticity using the re- 
sults of the metric functio ns re produced the Newtonian 
result derived by Ferra ro ( 1954 ) in the Newtonian limit 
(see Konno et al. 2000 ) justifies this formalism. 



A.3. Equations of order EqEb 

Finally, We consider the quantities ato and at 2 of order 
£b£q, which are related with the electric fields induced by 
stellar rotation. These functions obey the Maxwell equa- 
tion. Outside the star, we have 

_ x d 2 a m , 2e~ A da t 



dr 2 ' 


r dr 


h 2 2 


[(*■ 


f- - j w + J 


3 







dr 



^d a t2 2e A da 
dr 2 
2 

" ~3 



12 



e~ X a' 



3r 2 



uoa^, 











[(*- 






e 











+ ^2 Ua ^ 



(A.16) 



(A.17) 
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where the function u) of first order in eq was discussed in 
detail by Hartle (1967). These differential equations can 



also be solved analytically (see Muslimov & Tsygan 1986; 
Konno & Kojima 2000 for the detailed calculations). In- 
side the star, the assumption of a perfectly conducting 
interior yields the expressions (Bonazzola et al 
quet et al 



1995), 



1992: Boc- 



a*2 



2 

-ft, 



+ c 2 , 



2 



(A.18) 
(A.19) 



where c 2 is a constant fixed by the junction condition at 
the stellar surface. In this case, ato and Q>ti are not con- 
nected smoothly at the stellar surface, owing to the in- 
duced surface charge. 
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